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ABSTRACT 



Damage sustained by equipment used in a moving carrier 
suddenly changing its velocity can sometimes be attributed to in- 
adequate protective cushioning. Proper design of the protective 
cushioning can be realized by reducing the equipment and components 
to idealized mechanical systems and studying the responses as the 
system parameters are varied. This study considers the response to 
velocity shock with elastic impact. A set of design curves showing 
maximum acceleration of a flexible and vulnerable element of an 
equipment and the maximum deflection of the equipment's cushioning 
are produced using an analog computer. Additionally the computer's 
results are compared with analytical results for a few particular 
cases . 

The writer wishes to express his appreciation for the advice 
and assistance of Professor R. E. Newton in the selection and per- 
formance of this study . 
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SYMBOLS 



a k < 4 . • 

b Radian frequency of vibration of the packaged item n^on 
damped cushioning. 

C| Damping coefficient of linear cushioning of element rrij 
of packaged item. 

C 2 Damping coefficient of linear cushioning of packaged item 

• 

e Naperian base 2.718... 

f 

g Damped circular frequency of vibration of element mj of the 
packaged item. 

H, Constant defined by Eq. (28a) . 

Constant defined by Eq. (28b) . 

Constant defined by Eq. (23e) . 

k, Spring rate of linear elasticity of element m ( of packaged 
item . 

k 2 Spring rate of linear elasticity of packaged item m 2 . 

M Constant defined by Eq. (23h) . 

m, Lumped mass of element of packaged item. 
m f Lumped mass of packaged item. 
m 3 Lumped mass of outer container. 

N Constant defined by Eq. (23i) . 

P Constant defined by Eq.(23g) . 

Q Constant defined by Eq.(23f). 

t Time coordinate . 
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Time of rebound . 

Time coordinate after rebound . 

Initial velocity of m ( and m^ . 

Displacement of m , . 

Maximum acceleration of mj . 

Displacement of m A . Before rebound this is equal to isolator 
deflection . 

Displacement of linear cushioning of packaged item. 

Fraction of critical damping of an element m , of the packaged 
item . 

Fraction of critical damping of packaged item. 

Mass ratio yy\ t j yrv ; * 

frequency of m ( . 
frequency of . 



Undamped natural 
Undamped natural 
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I . Introduction 

One of the causes of mechanical damage in the transportation of 
packaged articles is inadequate protective cushioning. However , it is 
possible to predict what cushioning is required to protect the flexible 
and vulnerable elements of a given equipment. This prediction is based 
on the supposition that the container for the equipment, the protective 
cushioning, the equipment and the elements of the equipment can be 
idealized as a system of rigid masses, massless linear springs and 
linear viscous dampers (dashpots) . 

Figs. 1 and 2 illustrate the schematic representation of a package 

p 

and its representative mechanical system. 




Fig. 1. Schematic representation of a package. 
1 . Element of packaged item 
2a. The packaged item as a whole 
2b. Cushioning 
3. Outer container 
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Element of packaged Item 




\\\\\\\\\ \ \ \"\ V \ \ \ \ \ Support 
Fig. 2. Idealized mechanical system representing a package. 

In this study we shall be concerned with the responses of the 
package that is undergoing velocity shock with elastic impact. Velo- 
city shock is that shock which results from a sudden change in the linear 
velocity of the equipment or its support. If the outer container m., does 
not rebound from the support we call this behavior inelastic impact . On 
the other hand if mass m-j is negligibly small it will separate from the 
support when the net force applied by spring k , and dashpot c'n becomes 
tensile. We call this latter behavior elastic impact, or perfect rebound. 
The idealized system we shall use exhibits perfect rebound and is shown 
in Fig. 3. 
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Fig. 3. Idealized mechanical system consisting of a massive packaged 

item m 2 and a small element m^. The outer container is assumed 
massless . 

We now have a system with two degrees of freedom and , to make 
this system tractable, we will consider it to be linear. The free body 
diagrams of the two masses are shown in Fig. 4 for u = u = o. 





f ki ( x* - Xa ) 


k.(x,-x*) i 




( k*x* 






im*. 




w, 






. c* X* 









Fig. 4. Free body diagrams for and m^ of the idealized mechanical 
system. 

The differential equations of motion for m^ and m 2 are: 

< ^ 

w, X, = -k,(x,- X*) - c, ( X, - xj. (1) 

"W\ x Xa. = k,(X,'Xj + — Xa^~KaX z “C a X* . (2) 

It can be seen that there are six parameters: m^, m k^, kg, c'j and eg 



3 










. 



- ! i 






It is convenient to define 






L 
w> 



CO, 






J 



ft 



l\ 



, ) 



1 ’ z 
where the new parameters have the meanings: 

CO j , the natural undamped frequency of vibration of mass on spring k^ 
(x) 2 • the natural undamped frequency of vibration of mass mg on spring k , 
j , the fraction of critical damping for the subsystem consisting of 
m j , k j ; c ^ • 

g / the fraction of critical damping for the subsystem m g, k„, c^. 

Because both the length and time units may be chosen arbitrarily, 
it is possible to ignore absolute values of masses and frequen- 

cies CO i , U 2 . We use instead the mass ratio jk - m 1 /m„ and the fre~ 
quency ratio C0 S • The differential equations of motion become: 

X, = - wt (X,-Xj -a.fi.lo, (4) 



OO €P CO 

Xa. ~ “ Xx ” CUa 'X.a, ■“ X / . (5) 

For perfect rebound the first two terms on the right side of Eq. (S) do not 
apply after loss of contact with the support. Eq. (5) then becomes: 

L = -/6 X, C5“) 

From a practical standpoint is normally small and yw. is 
usually small. If we fix at a small value, say 0,01, and use the 
limiting form assumed by Eq. (5) as Jk approaches zero, we can reduce 
the number of system parameters to two: A and OJ t / *> x • 
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The behavior of the system thus idealized is studied in Section II and 



Appendix B . 

A second two-parameter idealization results from neglecting 
damping anc * confining attention to the effects of mass ratio 

and frequency ratio . This system is studied in Section II and 

Appendix B . 

With just two parameters, analytic solutions by known methods 
are feasible. However, such solutions are so laborious that it is im- 
practical to produce a set of design curves covering a range of values 
of these parameters. Accordingly, a set of design curves was obtained 
using an analog computer. 

CO 

The design results include the maximum acceleration X, 
of the element YM, of the- packaged item and the maximum deflection 
of the cushioning. The maximum acceleration of W\ ( is of interest 

cc 

because it is a direct measure of the maximum force Vv\ l X| fftAX acting 
on the structure of the element. The maximum deflection Xa-w\ax. is 
of interest because it determines the required physical dimensions of 
the cushioning. 
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II . Experimental Procedure 

The free body diagrams for vn and w\, of the idealized mechanical 
system shown in Fig. 4 of the Introduction are repeated below. 




C,( X x - X.) 



k,(x,-xj 

J. 




at X = 0 ; 

X, = X* = 0 3 

X, = x* = o. 



Fig. 4. Repeated. 



The differential equations of motion are Eqs. (4) and (5) of the 
Introduction . 

Xr -><=.). «> 

C C c cc 

X a - JX x, , (5) 

By setting: 



~ = X, - X 



( 6 ) 



~ n (i. U3, 



(7) 



i ; 

Eq. (4) becomes: 

~ J 1 ~ ~ J I \ VA ”T J-j’ j • 

Eqs. (5) and (7) are set up on the analog computer as indicated in 

Fig. 5. Equipment used is listed in Appendix E. 

The computer was scaled on a real time basis. Coefficient 

potentiometers were used to correct for the differences between nominal 
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condition 

h-cj 
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and actual values of resistors and capacitors. The Donner Comparator 



Bridge was used for precision measurement of the potentiometer settings. 
A diode limiter, type 6AL5 , was connected across the feedback of #1 
summing amplifier to cut off the output voltage when the voltage changed 
sign from minus to plus. The time of cutoff is the time of rebound. 

In making the "runs" the initial condition voltage was 

connected across the #3 operational amplifier. For the "runs" with 
damping and zero mass ratio, the coefficient potentiometer was set 

at 0.02, the coefficient potentiometer was set as required (0.0, 

0.1, 0.2, 0.4, 0.6 or 0.8) and the ju coefficient potentiometer was set 
at zero. For the "runs" with no damping and significant mass ratio, the 
and coefficient potentiometers were set at zero and the 

coefficient potentiometer was set as required (0.01, 0 . 1 , 0 . 2 , 0.4 or 
1.0). In ail "runs" the CJ Z coefficient potentiometers were kept at a 
constant value of 0.1 rad. /sec. and the Cd, coefficient potentiometers 
were varied to provide the desired frequency ratio ( CO, / C 0 Z ) . 



III. Analog Results for Velocity Shock with Perfect. Rebound 

The analog results are presented in the form of design curves , 
Figs. 6 to 13. The curves are plots of: 

(a) Xi ,wax/ Xo 00» / with frequency ratio as abscissa. 

X) max / & 0 GJ, is dimensionless and is called the '"shock 

O© 

transmissibility " . It is the ratio of maximum acceleration X t vr\«* 

o 

to the maximum acceleration X 0 which would result with no 

cushioning . 



(w X \ max / CO* , with frequency ratio as abscissa. 

CO 

Xi cnax / ^Co 6J* is called an "amplification factor". It is the 

ratio of the maximum acceleration of rrij to the maximum acceleration 
of m 2 with no damping in the cushioning. 

(c) X* / with frequency ratio as abscissa. 

Xx ma»y/( X /U> ) ls ca U ec * a "deflection factor". It is the ratio 

of the maximum deflection of the cushioning to the maximum deflection 



of the cushioning with no damping and negligible mass of the element m 



1 



(d) Xa. mayf^o/tOx.) with damping fraction f?> z as abscissa. For mas 
ratio = © , Xa max^/C k c / CO 3 .) is independent of the frequency 
ratio but is dependent on which can be seen from Eq„ (22). 

These curves can be utilized in the selection of adequate cushion- 
ing for a flexible and vulnerable element of equipment undergoing a 
velocity shock. 
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Fig. 7. Shock Transmissibility for linear undamped cushioning 
and significant mass ratios with perfect rebound . 
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Fig. 0. Shock. Transmissibility for linear damped cushioning 
with perfect rebound. Analog values are shown. 
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and significant mass ratios with perfect refco.ind 
Analog values are shown . 
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IV. Comparison of Analog and Analytical Result: 

The analytical results of Appendix C are cor p me . with mdoj 
results in Figs. 14 to 18. 

Figures 16 and 17 show the response c~r e ioi two i~n. 
made with the same values for the parameters. There is a diUerence ot 
approximately two per cent in the maximum magnitudes of the two runs 
In all cases the analog plot lags the analytical point plot The 
lag varies from plot to plot and ranges fiom one to three seconds m 30 
seconds . It is difficult to determine a precise comparison between mag- 
nitudes of the analog and analytical results,. This difficulty is *o 
the time lag and to the uncertainty of the exact analytical maximum mag- 
nitude. A comparison of the maximum magnitudes of the an dog results 
with the largest analytical values suggests that the error is within the 
range of plus or minus four per cent. 

A second method of estimating error is to use energy methods By 
equating the kinetic energy of the entire system before impact - o the 
potential energy stored in the cushioning at its maximum deflection the 



possible value can then be compared to the maximum value from the 



maximum possible value of 




- can be determined. This maximum 



analog results. With no damping 



( 




the comparison is as follows: 




£ 
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before 



We know that with yk * /. 0 , m,* kl x * = X x X a 

impact and CO* - V . 4 | . Substituting the ak.e in be energy 

Eq. (8) gives: 



± ^ 



a 3- 

X, 



o -3. 

k ^x 0 



= i to. 



o VL 



x 0 - 



^ a. A ' a (vv-ay 



4- cor va :l x= 



3L 



X S' 



Rearranging gives: 



£ - 



x 



X fyy^+M* 



Xo/ 



Cj 



A. 



and 



X 



•1 /YV-C^f 



Xo / to. 






which if 



X 



cl rvtv^4 



the maximum possible value for 7r ~77 r for this case. Wren we 

^ X. 

compare this value of JX to the maximum analog result we note 
that the analog result of 1.45 is about 2.6% high. The actual error is 
probably a little higher because some potential energy may be stored 
in the spring |( ( . 

A third method of estimating error is to compare the analytical 
and analog results for the maximum deflections of the cushioning with 
damping and zero mass ratio. The maximum deflection of the cushioning 
occurs before rebound and is dependent on 5^ as shown in Eq. (i,2). 
Comparison of the analytical and analog results shown in Fig. 18 
suggests of an accuracy of ± 2%. 
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Of further interest is a comparison with the tnalo j re^uk.'. pub 
lished in "Dynamics of Package Cushioning" fc/ R D, Mindlir who 
uses an amplification factor which is practically equal to tie impliiicii 
tion factor used in this study. 

For perfect rebound with: 



A = o-oi , 

Ai = - 

LL — 0.0 , 



Frequency ratio 

C a , j CO 


Amplificatior 

“^a,) rt\8^ 


i Factor 

oc 


’x t 

Dauber 


^ y , i -■*- q ^ 

Mindlin QQ 


1 


1.54 


1 . 56 


1.4 


1.65 


1.72 


2 


1.61 


1.69 


2.5 


1.5 


1.57 


3 


1.42 


1.45 


4 


1.24 


1.23 


5 


1.08 


1.07 


6 


1.15 


1.14 


7 


1.14 


1.13 


8 


1.12 


1 . 08 


9 


1.06 


o 

o 



The results compare favorably except at the low frequency ratios 



of 1.4 to 2.5 where the difference is approximately 4-5%. 
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Another comparison of Interest can be made with analog results for 
inelastic impact published in "Vibration and Shock Isolation JC by C. E. 
Crede [2]. Crede reduced data originally obtained by R. Do Mindlin |lj , 
At high frequency ratios the maximum acceleration of m occurs before re- 
bound so it is reasonable to expect the results for inelastic impact to be 
in close agreement with the results for perfect rebound. For ; 



= 0 . 0 / 



Frequency ratio 

to | / 




Tran suds si 
Crede 00 


bility 

Dauber 


3 


0.05 


0 „ 46 


0.46 


3 


0.1 


0.48 


0.46 


6 


0.05 


0.19 


0.178 


6 


0.1 


0.13 


0.18 



Crede' s values are taken from a plot and are not to be considered the 
exact experimental values . It can be seen that the agreement is reasonably 
close. 

The above comparisons of analytical and other analog results with 
the analog results of this study suggest that the accuracy of the latter is 
approximately t 4% . 
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Fig. 14. Comparison of analytical solutions with analog results. 
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Fig. 15. Comparison of analytical solutions with analog results. 
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Fig. 16. Comparison of analytical solutions with analog results. 
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Fig. 17. Comparison of analytical solutions with analog results 
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V. 



Conclusions 



Design curves have been obtained for piedic in ; ve ue.- - or o: 
a two degree- of -freedom linear system with perfect re bo n i T. » e?e 
curves are believed accurate within 1 The following ;ene. *1 feituret 

of the response behavior may be deduced 

a) Transmissibility curves show that in order to benefit significantly 
from cushioning the frequency ratio h J(J A must he greater than 

2 (Figs. 6 and 7). In this region increasing to, / Uj a decreases 
X, . The addition of small amounts of d a?rpln :j , | . f 

in the cushioning has negligible effect and greater damping Increases 

r 

X t (Fig. 6) . For a number of element * o: the itme ire- 

quency k , the acceleration a, { will be least for the most 
massive (largest value of m ) (Fig. 7). 

b) The maximum cushion deflection Xj , varies inversely with 

COjl and decreases monotonically with increasing damping 

(Fig. 12). 
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APPENDIX A 



Summary of Experimental Data 



The summary of experimental data is presented in Tables 1 to 

12 . 

Tables 1 to 6 list the maximum responses with damping and zero 
mass ratio . 

Tables 7 to 12 list the maximum responses with no damping and 
various mass ratios. Tables 11 and 12 are for the same mass ratio 
s 1.0 and they show the uncertain scatter established during 
separate runs . 
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Table 2 
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frequency 
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.08 


54 

66 


52 

52 




.1 


72 


52 




.12 


76 


52.5 


6 

n 

3 


.14 


77 


53 


.16 


76.5 


53.5 




.18 


76 


54 


O 

O 

J! 


.2 


73.5 


54.5 


.25 


67 


55 




.3 


61 


55.5 


d = 0 . c 


.35 


55 


56 


P-u 


.4 


50 


56 




.45 


54 


56 


H ^ t . V . 


.5 


54.5 


56 


} K 


.55 


54 


56 




.6 


53.5 


56 




.65 


52 


56 




.7 


50.5 


56 




.75 


50 


56 




.8 


51 


56 




.85 


51 


56 




.9 


51 


56 




.95 


50 


56 




1.0 


50 


56 
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Table 10 
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Table 11 





OJ, 


;,o 

max. 






.02 


16.2 


52 




.04 


30.5 


52 


c 


.06 


42.5 


53 




.08 


48 


54 




.1 


50 


57 




.12 


50 


63. a 




.14 


48 


65 




.16 


46 


68 




.18 


44 


70 


o.o 


.2 


42 


72 




.25 


40.5 


7*n 


(L s G.o 


.3 


42 


72 


P>- 


.35 


41 


72 




.4 


39.5 


72 


ik - l.o 


.45 


39 


72.5 


/ 


.5 


39.5 


72.5 




.55 


39 


72.5 




.6 


38 


73 | 


(first series 


.65 


38.5 


73 


of runs at 


.7 


38.5 


73 


fA =1.0) 


.75 


38 


73 




.8 


38 


73 




.85 


38 


73 




.9 


38 


73 




.95 


38 


73 




1.0 


38 


73 



Table 12 





GJ i 


CO 








X, 


X* 






.02 




14 


51 


t 


.04 




28.8 


51.5 


coX o ~ i:> " 0 


.06 




41 


52 




.08 




47 


54 


OlV- c, ‘ 1 


.1 




49 


56 


.12 




49 


GO 




.14 




47.5 


64 


£,= C '° 


. 16 




45.8 


67.5 


.18 




43.5 


69 . 5 




.2 




41.3 


70.5 


^ =0.0 


.25 




39 


71.5 


l 


.3 




41 


71.5 




.35 




40.6 


72 


o 

II 


.4 

.45 




39 

38.2 


72 

72 




.5 




39 


72 




.55 




38.8 


72 


(second series 


. 6 




38 


72 


of runs at 


.65 




37.8 


72 


/t =1.0) 


.7 




38 


72 


.75 




38 


72 




.8 




37.5 


72 




.85 




37.4 


72 




.9 




37.5 


72 




.95 




37.4 


72 




1.0 




37.3 


72 
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Appendix B 



Analytical Solution to the Differential Equations of 



Motion 



1 . Analytical solution for and no damping . 

With no damping Eqs. (4) and (5) become; 

X, = - U)/- ( X, - XJ . 



X x = X x - 



/<■ X. . 



The limiting form of Eq. (10) when approaches zero is; 

vo 

x„. = -oCX x . 

The initial conditions at t-o are; 

X, = X A = O ; 



X. - X, = X. 



Substituting these initial conditions in the general solutions of Eqs, 
(9) and (11) gives: 



X,= x. 



\ UV 



; t /l l ( l. ^ j 4 



ClJ \ 



co 



Y Ci ^ co, 



t 









AAav., U)*X , 

Eqs. (12)- (14) are valid up to the time of rebound of m«». Rebound 

occurs at time t , when the acceleration of m-y is zero. 

1 c 



41 



t 



IT 



l c- 



Xr° 



03 



(15) 






After rebound: Losing contact with the support requires that we 
replace Eq. (11) with : 

OC' 

X* = o (ir) 

The initial conditions are determined at the time of rebound, t,=Tr/^* 
We will start a new time coordinate with "t - O . The initial conditions 



are: 



X, = - ><c 

■c'-o 



CO. 



O), 3 -- 



OJ 



^JL /y^y 1 1"" 



03 1 
U, 



X 



« ,/ 
t - O 



- X / 1 + /C^kL- TT Xl ) , 



cO j uj ^ 



03, 



U, 



X 



<a- , . 
t- - 



o 



X 



a. 



t O 



= - x e . 



Substituting these initial conditions in the general solutions of Eqs. 
(9) and (11*) gives: 



*1 


> v 

= -^X c OA i . /Cfr^ ITU), ( 


> o)X+TTU),L ^ 


(16) 


*)<*), 






CO 

X| 


- <2TX C 03X 03, 

C'C-CC}' 


jCb*L TT 03, 




(17) 
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X, 



c« 

x, 



-t 



Eqs. (16) - (18) are valid for "t. o and 4 ~ ft/ „ 

2. Analytical solution for yc «i with damping. 

The differential equations of motion are equations (4) and 
of the Introduction . 

X. I — ~ X* i ( Aj ~ ) U) i C X ( “ X% ) • ( 2 



X. 



I, 



A. /X, O,' 



X. i 



ju ' '■*= 






X*. i 



(5) 



The limiting form of Eq„ (5) when yAt approaches Eero is ; 



= -<*. ; x 



a r-J c..'j X 



3L 



(IS) 



The initial conditions at "t-o are: 



X, = X A 



o 



= K = 



Substituting these initial conditions in the general solutions of Eqs, 

(4) and (19) gives: 

tv 



X , s e 4 (K u , ia, - |t - Q .x * ^ t j + Wt + $ „ 
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X, i ^ - v A jj ^ Q ) e ^ ^ c r ~t 

+ ( ” J* K„. - u.~y <;- C'i. -4- -§ ^ \\ u ) Ai.-k, g^ 

4 ( ~ -^cu h\ > 4 CX 3 ‘ C\ — Cx ) c* X'xa, b~f 



4* ( - \/~ i J -»• •* *-b ^ + a./* e 1 , <ct Lt . 


(21) 




- _X t e /u^ \y~t * 

b 


(22) 


where 


dx. — y~> a ^ ^ ^ x - 


(23a) 




b= i-frt) % . 


(23b) 




i = A &>, • 


(23c) 




r = <h ( /-a 3 ’)^ - 


(23d) 




K, = *. + ^<S. -ia -bP . 


(23e) 


<x=\ 


M (. “ V/ 4 v.r ~x<cf ■+*(*>,*') + |\) ( xoJb -xk> 1 ) 

i -■ b" ■» cv*- - 2. cc£ +• u). 3 ‘^ a 4- (. Xc_ b — x M\) 3 ' 


(23f) 


P - 


/V (— b"" >•(.</•— CJ, ) — - /V\ ( -2 fL,b “ -X S b* , 
(- PS u/'-aj-o. x- - ( «a.ft.b “ A^-b)^* 


(23g) 




(A = x, -f * 


(23h) 
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N = 



X. 

i 

b 



U l '''■c 



b 



Eqs. (20) - (22) are valid up to the time of rebound which occurs 
at time t., when the acceleration of mg is zero. 



b = j. [jl - ^ u f^.)] 
^■° b L ^ \ Jj 



(24) 



After rebound: Losing contact with the support requires that we 
replace Eq. (19) with: 



X*=° 

The initial conditions are determined at time of rebound 



i: 



’) 



We 



will start a new time coordinate with o , The initial conditions 
will be designated as follows: 



t'sc 



e 

X 



1 



“ *»*, , 


' -fc 2=0 


= *■*, , 


M -t*o 



■> 



"X: 



Substituting these initial conditions in the general solution of Eqs 
(4) and (19‘) gives: 

Xt— cr, ^ ( H( < L + H x ij- X sl.^ t 



25 ) 



^■i - ( “i hi ) e ^ ixj ^_ 

+ ( ihl, - 9 ‘H, 

X..U = 



X:-. f X*. t . 

•L. -r 



( 27 ) 
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where 



Hi ~ X 't, “ Xa -t, - 

V 

Eqs. (25) - (27) are valid for "t^o and . 

3. Analytical solution for _ D y jz o 

With no damping Eqs. (4) and (5) become: 

Xi =; — i C X i — Xx V 

'v/ SL. . u»t 

— CO^. Xjl X * . 

The initial conditions at "t = o are: 



|Z«aji 



(28b) 



(10) 



x v = X 



a. - ° j 



- X*. = x. 



(a) Because the solution of Eqs. (9) and (10) becomes intractable in 
literal form, the following numerical values are assigned: 



Co ( = 0 . x 




><. = S-o o 


GJ ^ --0. 1 




yU_ - O. X . 



Substituting the initial conditions and the above assigned values in the 
general solutions of equations (9) and (10) gives: 

X, = - . a-iilir •+• C.C.STS-, x, . (29) 



X, — h> - s* 3 . x . OW44 1* 



(30) 



46 



= job.sr t^.x-xbieT + Xhisy^ „ ozi44 t . 

Equations (29) - (31) are valid up to time of rebound of . Ke 

bound occurs at time "t, when X z is zero. Time T i is determined by 

a trial and error procedure and 

"fc, — 35.35" seconds. (32) 

1 **.-c 

After perfect rebound: Losing contact with the support requires 
that we replace Eq. (10) with: 

SU = /•.u,* ( x, - X.*.) . ’ do') 

The initial conditions are determined at time of rebound t. t . We will 
start a new time coordinate with t'^o . The initial conditions are: 





= -STS-S , 










'*'*• G's. c> 


- 



=. -^7.3 . 

Substituting these initial conditions in the general solutions of Eqs . 
(9) and (10’) gives: 

X,r + 6 yu+.a/f/t'. (33) 



oo 

X, - 



SLO.C. . SU/ f / t ' + A.O. /Zixu . SL/f/ 7?-'. 



(34) 



X«a. - -?3>2L ft 4:' (35) 

Eqs. (33) - (35) are valid for and 1 i — 35.35 seconds. 

(b) For additional comparison with analog results the mass ratio yx, 
is set at unity and the assigned values for oj ( ^ and X a 

remain 0.2, 0.1 and 500 respectively. Substituting the initial conditions 
into the general solutions to (9) and (10) gives: 
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X, = . 5 , A . . ?.i *.* » 71 * l. TVj/., .Mr- t (36) 

X, “ .o~< • , JJ.V'// - -• 3£, 2 , , .CH-t'iZ (27) 

« //lit* . - v i / "/■' •/- .,aC v <c6rvT\- (38) 

Eqs. (36 - (38) are valid up to time of rebound which occurs at time t ( 
when is zero. 

*ti = *f-C.o or seconds. (39) 

^ u 

After perfect rebound: Losing contact with the support requires 
that we replace Eq. (10) with Eq. (10'). The initial conditions are 
determined at time of rebound ~i ' , . We will start a new time 
coordinate with "t - o . The initial conditions are: 






X| -e'=o 






o 



O 



X, 









Substituting the initial conditions into the general solutions of Eqs. 
(9) and (10') gives: 

X, = -^ 7.7 ---Hi,/ 1' -tfi.t *lU . cu'-x * *" -*■?.*#*« (40) 

X, “ 7 c*«f . . MiTjLef T' ; (41) 

~ <77.7 - ‘7"7J-«’X’ ’£ ,.ZX:* i 2T ir ' ^ ^4 2) 

Eqs. (40) - (42) are valid for "t-'s-c and t is 45.05 second". 
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Appendix C 



Summary of Analytical Data 



Table 13 






Analytical data for: 






j2>, = o.o/ 


CUJ / = 


o.z. 


& - O . u L 


OUa. — 


o./ 


£L -O.o 


>< = 
e> 


-6"o o 






The solutions are ( from Eqs . (20) - (22) and (25) - (27^ ) 
Before x'ebound 

y* i „.oou.t , i 

X, - 33.1 /Cfwu .JX't c± 



. . , O j. 1* i /i — .o«^.b» . -J. 

33 i <c; ® 



V . - <r io~l ,oq? ~t 



After rebound 



ct 



X , - H.X <X 






,C&ZJ . 



-h 



a -oo-x^-t t O-~t 



t: , 
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Table 13 
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Table 14 



Analytical data for: 



X =■ <•*<' 






X <• 



- O. X. 



fX, <>• > 

<u ^ ^ / 

^C.. «. 6'Of 



The solutions are Eqs . (30), (31), (34) and (35): 



Before rebound 

/Ci -■ -u.x „?a- . ■-< Wi' - - * . x .otfi'-H'i -■ ^ 

X,a.~ X. ...... > . i St- t -f -S>- /5xV/ dd' t /4 J f’£' * (31) 



After rebound 



.Go 

Xi - ^ . . . 


. /•// r ' v ^ v... -/ ..c^ xiqjt , 


(34) 


X-, > i. 


^ /-b'C ^ 


- 




f j.'l, c c-^ .Wf/ ^ } 


(35) 
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Table 14 



t 


t' 


cc 

X, 


t^X*. 


0 




0 


0 


4 




- 2.0 


19.5 


8 




-14.1 


36.0 


12 




-37.4 


47.3 


16 




-61.7 


52.3 


20 




-72.7 


51.0 


24 




-61.6 


43.8 


28 




-32.3 


31.7 


32 




+ 1.6 


15.5 


35.35 


0 


+ 22.4 


- .5 




2 


+ 29.0 


- 9.7 




4 


+ 30.2 


-19.5 




6 


+ 25.6 


-29.6 




8 


+ 16.3 


-40.0 




10 


+ 3.8 


-50.4 




12 


- 9.4 


-60.9 
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Table 15 



Analytical data for: 



II 

X 


<9, <9 


CU, - 


O.JL- 


ll 




SL. ~~ 


O 4 / 




/^O 


r 

Xo = 


<S"00 



The solutions are Eqs . (37), (38), (41) and (42): 

Before rebound 

X \ ~ Z? • £ 34, «X x2 S” 'i' 4 (37) 

X-J. = /Xx. iT -9-Xx/ £ (38) 

After rebound 

'Xi = oJU^ . ^7.0 /t&L* . xzxz ~t , (4^ 

,/ 

X a. = -Z7,7 - X% . 7 £' -r /1 9. ! ai^.xsxz ~fc + ?77. 7/^ ,'jzxf£>> (42) 
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Table 15 



t 


1 

t 


X, 1 


J 

u / 

• 41 I 


0 




0 


0 


4 




- 2.0 


19 „ 5 


8 




- 12.8 


36.3 


12 




- 29.6 


49.6 


16 




- 40.7 


59.5 


20 




- 39.1 


66 . 3 


24 




- 30.5 


68 . 8 


28 




- 26.0 


65.2 


32 




- 28.8 


55.5 


36 




- 30.2 


41.6 


40 




- 20.7 


25.8 


44 




- 2.2 


9.0 


46.05 


0 


+ 7.0 


- .6 




2 


+ 13.6 


- 9.1 




4 


+ 15.9 


- 18.8 




6 


+ 13.3 


- 29.0 




8 


+ 6.6 


- 39.8 




10 


- 2.3 


- 50.8 
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Appendix D Operating Routine 

Preliminary preparations: The coefficient potentiometers were 
calibrated with the comparator bridge. This calibration was verified 
several times during the course of the study. The recorder's timer was 
calibrated with an electric stop watch. 

Daily preparations: The computer’s amplifiers were balanced 
after an initial warmup period of at least twenty minutes and at the 
beginning of each series of runs. The initial conditions were verified 
after each series of runs . The recorder channels were calibrated before 
each series of runs. The calibration was checked after the series of 
runs and recalibrated if necessary. Capacitors and resistors were cali- 
brated with a capacitance bridge and a standard resistor respectively. 
Each capacitor was inserted in a separate integrating circuit and checked 
for improper period and excess leakage. Coefficient potentiometers in 
the problem circuit were adjusted for any day to day variation of the 
capacitors ." A check off sheet was maintained for all daily preparations . 

In all runs the potentiometers were kept at a - constant value 

of 0 . 1 rad. / sec. and the potentiometers were varied to give the 
desired OJ, / uJ, ratio. The ^ ^ and jJ. potentiometers were set 
with their required values and remained constant for the' complete range 
of the OJj / U* ratio. The output of Amplifier #1 was always fed to 
one of the recorder channels to verify proper operation of the diode 
limiter. At cutoff by the diode limiter the output voltage of Amplifier #1 
was approximately 0.4 volts. The maximum output of Amplifiers #5 and 
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#6 were recorded on the data sheets „ A rough plot of X, and * 
versus was made concurrent with computer operation. This 

rough plot served as an indication of unexpected responses. 



ij 
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Appendix E 



Equipment 



a. Analog Computer, Model 30 
Donner Scientific Company 
Berkeley, California 

b. Comparator Bridge, Model 50 
Donner Scientific Company 
Berkeley, California 

c. Four Channel Recorder, Model 154-1100A 
Sanborn Company 

Cambridge, Mass. 

d. Vacuum Tube Voltmeter, Model 311 
Simpson Electric Company 
Chicago Illinois 
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